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Abstract 

We consider a family of pseudo differential operators {A + a a A a ^ 2 ; a <G (0, 1]} on W l for 
every d > 1 that evolves continuously from A to A + A"/ 2 , where a 6 (0, 2). It gives rise to 
a family of Levy processes {X a , a £ (0, 1]} in M. d , where X a is the sum of a Brownian motion 
and an independent symmetric a-stable process with weight a. We establish sharp two-sided 
estimates for the heat kernel of A + a a A a / 2 with zero exterior condition in a family of open 
subsets, including bounded C 1,1 (possibly disconnected) open sets. This heat kernel is also the 
transition density of the sum of a Brownian motion and an independent symmetric a-stable 
process with weight a in such open sets. 

Our result is the first sharp two-sided estimates for the transition density of a Markov process 
with both diffusion and jump components in open sets. Moreover, our result is uniform in a in 
the sense that the constants in the estimates are independent of a £ (0, 1] so that it recovers 
the Dirichlet heat kernel estimates for Brownian motion by taking a — > 0. Integrating the heat 
kernel estimates in time t, we recover the two-sided sharp uniform Green function estimates 
of X a in bounded C 1,1 open sets in R d , which were recently established in Q3] by using a 
completely different approach. 
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1 Introduction 

Many physical and economic systems should be and in fact have been successfully modeled by dis- 
continuous Markov processes; see for example, [27, 26, 28J and the references therein. Discontinuous 
Markov processes are also very important from a theoretical point of view, since they contain stable 
processes, relativistic stable processes and jump diffusions as special cases. Due to their importance 
both in theory and in applications, discontinuous Markov processes have been receiving intensive 
study in recent years. 
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In general, Markov processes may have both diffusion and jump components. A Markov process 
having continuous sample paths is called a diffusion. Diffusion processes in M. d and second order 
elliptic differential operators on M. d are closely related in the following sense. For a large class of 
second order elliptic differential operators C on M. d , there is a diffusion process X in R rf associated 
with it so that C is the infinitesimal generator of X, and vice versa. The connection between C and 
X can also be seen as follows. The fundamental solution p(t, x, y) of dtu = Cu (also called the heat 
kernel of C) is the transition density of X. Thus obtaining sharp two-sided estimates for p(t, x, y) 
is a fundamental problem in both analysis and probability theory. In fact, two-sided heat kernel 
estimates for diffusions in M, d have a long history and many beautiful results have been established. 
See |21t I23j and the references therein. But, due to the complication near the boundary, two-sided 
estimates on the transition density of killed diffusions in a domain D (equivalently, the Dirichlet 
heat kernel) have been established only recently. See [22\ [23| [2^] for upper bound estimates and 
|31j for the lower bound estimate of the Dirichlet heat kernels in bounded C 1,1 domains. 

The infinitesimal generator of a discontinuous Markov process in W 1 is no longer a differential 
operator but rather a non-local (or integro-differential) operator C For instance, the infinitesimal 
generator of a rotationally symmetric a-stable process in M. d with a S (0, 2) is a fractional Laplacian 
operator A"/ 2 := —(—A)"/ 2 . Most of the recent studies concentrate on pure jump Markov pro- 
cesses, like the rotationally symmetric a-stable processes, that do not have a diffusion component. 
For a summary of some of these recent results from the probability literature, one can see [H [7] 
and the references therein. We refer the readers to [U 16] for a sample of recent progresses in the 
PDE literature. 

Recently in [9], we obtained sharp two-sided estimates for the heat kernel of the fractional 
Laplacian A a / 2 in D with zero exterior condition (or equivalently, the transition density function 
of the symmetric a-stable process killed upon exiting D) for any C 1,1 open set D CW 1 with d > 1. 
As far as we know, this was the first time sharp two-sided estimates were established for Dirichlet 
heat kernels of non-local operators. Since then, studies on this topic have been growing rapidly. In 
PHI [HI [12] , the ideas of [9] were adapted to establish two-sided heat kernel estimates of other pure 
jump Markov processes in open subsets of M. d . In [19], the large time behaviors of heat kernels for 
symmetric a-stable processes and censored stable processes in unbounded open sets were studied. 
Very recently in [21 [3], the heat kernel of the fractional Laplacian in non-smooth open set was 
discussed. We refer the readers to [8] for a survey on the recent progresses in the heat kernel 
estimates of jump Markov processes. 

However until now, two-sided heat kernel estimates of Markov processes with both diffusion 
and jump components in proper open subsets of R rf have not been studied. The fact that such a 
process X has both diffusion and jump components is the source of many difficulties. The main 
difficulty stems from the fact that such a process X runs on two different scales: on the small scale 
the diffusion part dominates, while on the large scale the jumps take over. Another difficulty is 
encountered at the exit of X from an open set: for diffusions, the exit is through the boundary, 
while for the pure jump processes, typically the exit happens by jumping across the boundary. For 
a process X that has both diffusion and jump components, both cases will occur, which makes the 
process X much more difficult to study. 

In this paper, we consider Levy processes that are independent sums of Brownian motions and 
(rotationally) symmetric stable processes in R rf with d > 1. We establish two-sided heat kernel 
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estimates for such Levy processes killed upon exiting a C ' open set. The processes studied in this 
paper serve as a test case for more general processes with both diffusion and jump components, just 
like Brownian motions do for more general diffusions. We hope that our study will help to shed new 
light on the understanding of the heat kernel behavior of more general Markov processes. Although 
two-sided heat kernel estimates for Markov processes with both diffusion and jump components 
in M. d have been studied recently in [17^ I29j . as far as we know, this is the first time that sharp 
two-sided estimates on the Dirichlet heat kernels for Markov processes with both diffusion and 
jump components in proper open subsets are established. 

Let us now describe the main result of this paper and at the same time fix the notations. 
Throughout this paper, we assume that d > 1 is an integer and a E (0,2). Let X° = (X®, t > 0) 



be a Brownian motion in M. d with generator A = ^2 



d 

=1 dx: 



.. Let Y = (Y t , t > 0) be a (rotationally) 



symmetric a-stable process in Mr, that is, a Levy process such that 



IE, 



for every x G 



l d and f G M d . 



The infinitesimal generator of a symmetric a-stable process Y in M d is the fractional Laplacian 
A a / 2 , which is a prototype of non-local operators. The fractional Laplacian can be written in the 
form 

A a/2 «W = lim / (u(y) - <x)) dy, (1.1) 

<40 J{ y( zRd. \ y - x \ >E } \x - y\ a+a 

where A(d,a) : = a2 a - 1 vr- d / 2 r(^)r(l - § )~ 1 . Here T is the Gamma function defined by T(A) : = 
J* °° t ■ 1 e — *d£ for every A > 0. Assume that X° and Y are independent. For any a > 0, we define 
X a by Xf := X® + aYf We will call the process X a the independent sum of the Brownian motion 
X° and the symmetric a-stable process Y with weight a > 0. The infinitesimal generator of X a is 
A + a a A a / 2 and 



E, 



for every x E R d and £ G M. d . 



Since 



a a A(d,a) 

cos ^ • y)) \ y \d+ a d y> 



the density of the Levy measure of X a with respect to the Lebesgue measure on 1R is 

J a (x, y) = f(\x - y\) := a a A(d, a)\x - y\^ d+a l 



The function J a (x,y) determines a Levy system for X a , which describes the jumps of the process 
X a : for any non-negative measurable function / on IR + x M. d x M. d with f(s, y, y) = for all y G M. d , 
any stopping time T (with respect to the filtration of X a ) and any x G M. d , 



Y,f(s,X-_,X-) 



[ T ( [ f( S ,X*,y)J a (X*,y)dy)ds 
J \JR d / 



(1.2) 



(see, for example, \15\ Proof of Lemma 4.7] and [16\ Appendix A]). Let p a (t, x, y) be the transition 
density of the process X a with respect to the Lebesgue measure on W d , which is known to exist 
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and is jointly continuous on (0, oo) x R d x R d . For any A > 0, the process (XX^_ 2t ,t > 0) has the 
same distribution as (ATf A(a 2)/a ,t > 0) (see the second paragraph of Section [21), so we have 

p^ a - 2)/a (t :X ,y) = \- d p a (X-H,\- 1 x,X- 1 y) for t > and x, y G R d . (1.3) 

The following sharp two-sided estimates on p a (t,x,y) follows from (jl.3p and the main results 
in [T71 [29] that give the sharp estimates on p 1 ^, x, y). 

Theorem 1.1 There are constants Ci > 1, % = 1,2, such that, for all a G [0, oo) and (t,x,y) G 
(0, oo] x R d x R d 

a a t 



y\d+a 



Cf 1 (V d/2 A (a Q t)" d/Q ) A ^f-^e-Cils-tfl 3 /* + ( a a t)~ d / a A 
< p a (t,x,y) < d (t- d ' 2 A {a a t)- d l°^ A ^-^g-lx-ylVCa* + ( a «f)-«*/« A . 

Here for a, 6 G R, a Ab := min{a, 6} and a V b := max{a, b}. In particular, we have 

Corollary 1.2 For any M > and T > 0, i/iere is a constant C3 > 1 depending only on M and 
T such that, for all a G (0, M] and (t, x, y) G (0, T) x R d x R d 

C7 l ( t -d/2 e -C 2 \ x - y \yt + t -d/2 A 



3 \ 1 \ x _ y \d+a 



< p a (t,x,y) < C 3 (t-d/2 e -\x-y\ 2 /C2t + t -d/2 A 
where C2 > 1 is i/ie constant in Theorem \l.l[ 



a a t 



\x - y\ d+a 



Recall that an open set D in M d (when d > 2) is said to be C 1,1 if there exist a localization radius 
R > and a constant A > such that for every z G 91?, there is a C ^-function 4> = <j) z : R^ 1 — > R 
satisfying 0(0) = 0, V0(O) = (0, ...,0), ||V0[|oo < A, |V^(x) - V<f){z)\ < A\x - z\, and an 
orthonormal coordinate system CS Z : y = (yi, ■ ■ ■ ,yd-i,yd) '■= (y, yd) with origin at z such that 
B(z,R) n D = {y = (y,y d ) G B(0,R) in CS Z : y d > <p(y)}. The pair (R, A) will be called the C 1 ' 1 
characteristics of the open set D. By a C 1 ' open set in R we mean an open set which can be 
written as the union of disjoint intervals so that the minimum of the lengths of all these intervals 
is positive and the minimum of the distances between these intervals is positive. Note that a C ' 
open set can be unbounded and disconnected, and that a bounded C ' open set have only finitely 
many connected components. 

For an open set D cR d and x G D, we will use 5d(x) to denote the Euclidean distance between 
x and D c . For an open set D C R d and (ro,Ao) G (0,oo) x [l,oo), we say the path distance 
in each connected component of D is comparable to the Euclidean distance with characteristics 
(ro,Ao) if the following holds for any r G (0, ro): for every x,y in the same component of D with 
Sd(x) A So(y) > r, there is a rectifiable curve I in D connecting x to y so that the length of / is no 
larger than Xq\x — y\. Clearly, such a property holds for all bounded C 1 ' open sets, C > open sets 
with compact complements and domains above graphs of C ' functions. 

For any open subset D C R d , we use to denote the first time the process X a exits D. We 
define the process X a ' D by X^° = XI for t < t% and X a t D = d for t > r% , where d is a cemetery 
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point. X a ' is called the subprocess of X a killed upon exiting D. The infinitesimal generator 
of X a ' D is (A + a a A a / 2 )\£). It follows from [17] that X a,D has a continuous transition density 
p a D {t,x,y) with respect to the Lebesgue measure. 

The goal of this paper is to get the following sharp two-sided estimates on p a D {t,x,y) for any 
C ' open set D in which the path distance in each connected component of D is comparable to 
the Euclidean distance. 

Let 

(l A ^M) (l A ^M) (t-^ e - c \*-y\ 2 /t + A r </*)) 

when x, y are in the same component of D, 

when x, y are in different components of D. 

One can easily show that, when D is bounded, the operator —(A + a a A a / 2 )\i) has discrete 
spectrum (see, for instance, the first paragraph of the proof of Theorem 11.31 (ii) and (iii) in Section 
4). In this case, we use \1' D > to denote the smallest eigenvalue of — (A + a a A a ^ 2 )\o- Denote by 
D(x) the connected component of D that contains x and let \^' D ( X ' > be the smallest eigenvalue 
of-(A + a"A a / 2 )| D(a;) . 

Theorem 1.3 Let d>l. Suppose that D is a C 1,1 open set in M. d with characteristic (R,A) such 
that the path distance in each connected component of D is comparable to the Euclidean distance 
with characteristics (ro,Ao). 

(i) For every M > and T > 0, there are positive constants Ci = Ci(R, A, ro, Ao, M,a,T) > 1, 

i = 4, 5 ; such that for all a G (0, M] and (t, x, y) G (0, T) x D x D, 

C^ 1 h a C5 (t,x,y) < p a D (t,x,y) < C 4 h a 1/C5 (t,x,y). 

(ii) Suppose in addition that D is bounded and connected. For every M > and T > 0, there is a 

constant C 6 = C e (D, M, a, T) > 1 so that for all a £ (0, M\ and (t, x, y) G [T, oo) x D x D, 

C^e- tx " D 5 D (x)5 D (y) < p a D (t,x,y) < C 6 e~* A ?'° 5 D {x) 5 D (y). 

(iii) Suppose that D is bounded but disconnected. Then for every M > and T > 0, t/iere are 
constants Ci = Ci(D, M, a,T) > 1. z = 7, 8, such that for all a G (0, M], t G [T, oo), the 
following hold. 

(a) Ifx,y are in the same component D(x) of D, 

C^e^^ 5 D (x)6 D (y) < p a D (t,x,y) < C 7 (eT* ^ + (l A (a a t)) e "* A " D ) fe(x) «J D (y). 

(b) Ifx,y are in different components of D, 
C 8 1 a a te- t ^' D(X)vX " D(y) n D (x)S D (y) < p a D (t,x,y) < C 8 (1 A (a a t)) e~* x *'° S D (x) 6 D (y). 



hc(t,x,y) := < 
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Remark 1.4 (i) Unlike the Brownian motion case, even though D may be disconnected, the 
process X a,D is always irreducible when a > because X a,D can jump from one component 
of D to another. When a > is smaller, the connection between different components of D 
by X a becomes weaker. The estimates given in Theorem 11.31 present a precise quantitative 
description of such a phenomenon. Letting a — >■ 0, Theorem 11.31 recovers the Dirichlet heat 
kernel estimates for Brownian motion in D (even when D is disconnected) ; see [20} [3T] and 
the reference therein for the latter. In particular, for x and y in different components of D, 
we have lim. a -+o+ Pi)(t,x, y) = for all x,y > 0, which is the case for Brownian motion. 

(ii) In fact, the estimates in Theorem 11.3( 1) will be established under a weaker assumption on 
D: the lower bounded estimate is proved under the uniform interior ball condition and the 
condition that the path distance in each connected component of D is comparable to the 
Euclidean distance (see Theorem I2.4p . while the upper bound estimate is proved under a 
weaker version of the uniform exterior ball condition (see Theorem I3.9[) . Here an open set 
D C M. d is said to satisfy the uniform interior ball condition with radius R\ > if for every 
x € D with 5d{x) < Ri, there is z x S 3D so that \x — z x \ = 5d(x) and B(xo,Ri) C D for 
xo := z x + R\{x — z x ) /\x — z x \. We say D satisfies a weaker version of the uniform exterior 
ball condition with radius R\ > if for every z G 3D, there is a ball B z of radius R\ such 
that B z C (D) c and 3B Z C\3D = {z}. 

Integrating the heat kernel estimates in Theorem 11.31 over time t yields the following two-sided 
sharp estimates of the Green function of X a in bounded C jl open sets, which were first obtained 
in [13] by a different method. We will not give the details in this paper on how these estimates 
can be obtained by integrating the estimates in Theorem 1 1.31 Interested readers are referred to the 
proof of [SI Corollary 1.2], where the sharp estimates for the Green functions of symmetric stable 
processes in bounded C ' open sets are obtained from the sharp heat kernel estimates for the heat 
kernels by integration over time t. 

Define for d > 3 and a > 0, 

{ \ x -y\d-2 ( 1 A "^""jr^J when x, y are in the same component of D, 

\ x -y\d-2 ( 1 A &D< ^y^ J when x, y are in different components of D; 

for d = 2 and a > 0, 

{log ^1 + SD \^-y\± V ' > ) when x, y are in the same component of D, 

a a log ^1 + -^j^Zyjr^J when x, y are in different components of D; 
and for d = 1 and a > 0, 

{((^D^cWy)) 1 / 2 A Sd ^^ d ^ when x,y are in the same component of D, 

° Q ( (°~D(x)°~D(y)) 1 ^ 2 A ^° jx-^j ) when x,y are in different components of D. 

Corollary 1.5 Lei M > 0. Suppose that D is a bounded C 1,1 open sei in M. d . There exists 
Cq = Cg(D, M, a) > 1 such that for all x,y E D and a// a € (0, M] 

Cg- 1 ^,?,) < G%(x,y) < C 9 g a D (x,y). 
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This paper is a natural continuation of [9], where sharp two-sided heat kernel estimates for 
symmetric a-stable processes in C 1,1 open sets are first derived, as well as [T3], where the boundary 
Harnack principle for X a is established. Some ideas of the approach in this paper can be traced 
back to [9] but a number of new ideas are needed to handle the combined effects of Brownian 
motion and discontinuous stable process. A comparison with subordination of killed Brownian 
motion is used for the lower bound short time heat kernel estimates for X a,D . We would like to 
point out that, unlike pj, the boundary Harnack principle for X a is not used directly in this paper. 
Instead we use one of the key lemmas established in [13] to obtain the upper bound of the heat 
kernel (see Lemma 13. X j) . Theorem ll.3f i) will be established through Theorem 12.41 and Theorem 
13.91 which give the lower bound and upper bound estimates, respectively. In contrast to that in 
[TU1 [TT1 112j . the proof of large time heat kernel estimates in Theorem ll.3( ii)-(iii) does not use 
intrinsic ultracontractivity of X a,D . The proof presented here is more direct, and uses only the 
continuity of A" ,£> and its corresponding first eigenfunction in a E (0, M], which is established in 
|18j . Lastly, we point out that the approach of [3] relies critically on the fact the symmetric stable 
processes do not have diffusion component and so it is not directly applicable to the processes 
considered in this paper. 

We will use capital letters C\, C2, . . . to denote constants in the statements of results, and their 
labeling will be fixed. The lower case constants ci,C2 3 ... will denote generic constants used in 
proofs, whose exact values are not important and can change from one appearance to another. The 
labeling of the lower case constants starts anew in each proof. The dependence of the constant c 
on the dimension d will not be mentioned explicitly. We will use ":=" to denote a definition, which 
is read as "is defined to be". We will use d to denote a cemetery point and for every function /, 
we extend its definition to d by setting f(d) = 0. We will use dx to denote the Lebesgue measure 
in W 1 . The Lebesgue measure of a Borel set A C M. d will be denoted by \A\. 



2 Lower bound estimate 



In this section, we assume that D is an open set in R satisfying the uniform interior ball condition 
with radius Ri > and that the path distance in each connected component of D is comparable to 
the Euclidean distance with characteristics (ro,Ao). Observe that under the uniform interior ball 
condition, the condition that the path distance for each connected component in D is comparable 
to the Euclidean distance is equivalent to the following: there exist r<i , A > such that for all 
r G (0, T2\ and all x,y in the same connected component of D with 5d(x) A 5d(u) > r, there 
is a rectifiable curve I in D connecting x to y so that the length of I is no larger than \\x — y\ 
and o~d{ z ) > r fc> r every z G I. The latter is also equivalent to the following, which is called the 
connected ball condition in [20]: For all r G (0, 7*2] and x,y in the same connected component of 
D with 5d(x) A 5d(u) > r, there exist m and x^, k = 1,2, ...,m such that xq = x, x m = y, 
Xfc-l G B(xk, §) C B(xk, r) C D and r ■ m < Xq\x — y\. 
Observe for all A, a > and £, x G W 1 , 



E, 



o *.(A(*« ■ -Xg)) 



o i(aXt)-(Y t/x2 -Y ) 



-t(|£| 2 +(aA( a - 2 )/ Q ) a |£| a ) 



It follows that if {X^ ,D , t > 0} is the subprocess in D of the independent sum of a Brownian motion 
and a symmetric a-stable process on M. d with weight a, then {AX°L^ t ,t > 0} is the subprocess in 
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XD of the independent sum of a Brownian motion and a symmetric a-stable process on M. d with 
weight a\( a ~ 2 ^ a . So for any A > 0, we have 

P\D a ~ 2)/a (t,x,y) = \- d p a D (\~ 2 t, X^x, X^y) for t > and x, y G XD. (2.1) 

The above scaling property of X a will be used throughout this paper. For t > 0, we define 

at:=at (2-»)/(2a)_ (22) 

This notation will be used in this paper when we scale an open D by s~ 1//2 to s _1 / 2 D. 

We first recall the definition of subordinate killed Brownian motion: Assume that U is an 
open subset in M. d and Tf is an a/2-stable subordinator independent of the killed Brownian motion 
X°' U . For each a > 0, let T a be the subordinator defined by T 4 a := t + a 2 Tf. Then the process 
{Z®' : i > 0} defined by Z"' = X^ 7 is called a subordinate killed Brownian motion in U. Let 
qfj(t,x,y) be the transition density of Z a,u . Then it follows from [301 Proposition 3.1] that 

pfr(t, z, w) > q?j(t, z, w), (t, z, w) G (0, oo) X U X U. (2.3) 

We will use this fact in the next result. 

Lemma 2.1 Suppose that M and T are positive constants. Then there exist positive constants 
C% = Ci(Ri,ro, Ao, a, T, M), i = 10, 11, such that for all a £ (0, M], t G (0, T] and x, y in the same 
connected component of D, 



p a D (t,x,y)>C w f 



- d/2 (lA S J^f) (lA^M) e -C^-v?/t, 



Proof. Suppose that x and y are in the same component, say U, of D. Let pu(t,x,y) be the 
transition density of the killed Brownian motion in U. It follows from [201 Theorem 3.3] (see also |31|, 
Theorem 1.2]) that there exist positive constants c\ = cx(i?i,ro, Xo,a,T) and C2 = C2(i?i, ro, Aq, a) 
such that for any (s,x,y) G (0, 2T] xU x U, 

PU(S, x, y) > Cl (l A Mfl) (l A ^M) ,^MV. 

(Although not explicitly mentioned in [20], a careful exam of the proofs in [20] reveals that the 
constants c\ and C2 in above lower bound estimate can be chosen to depend only on (R\, ro, Ao, a, T) 
and (i?i, ro, Ao, a), respectively.) Since p t -i/2jj(u, i~ 1//2 x, t~ l l 2 y) = t d / 2 pu(ut,x,y), we have for 
t < T and (u, x, y) G (0, 2] x U x U, 

p^ /2u ( Uj t-^ 2 x,t- l / 2 y) > Cl (l A f 1 A u-^eT^-^/M. (2.4) 

V ytu J \ \Jtu ) 

Let |U a *(n, s) be the density of a 2 T u , where at is defined in (|2.2p . Then it follows from the 
definition of the subordinate killed Brownian motion (for example, see [U page 149]) that for every 
1/3 < b < 1 and < t < T, 

POO 

tgLi/tufrt-WxJ-Wy) = / p^uis^-^x^-^yMb + a^neds) 

Jb 
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b 



oo 

p t -i/*u(s, t~ 1/2 x, t- 1/2 y)fi at (6, s - b)ds 



OC 







P t -V2jj{s + b, t~ 1/2 x, t- 1/2 y)/j at (6, s)ds. 
Consequently, by Q and (H3|), for every 1/3 < b < 1 and < t < T, 



> 



P a t t 1/2u {b,t- l / 2 x,r l l 2 y) 



> q% /2TI {b^l 2 x^/ 2 y) 
> 



t-l/2[/* 
1 

P t -i/2jj{s + b, t~ 1/2 x, t- l/2 y)fi at (b, s)ds 



> 



ci 
2 



> C3 1 



„ (lA 6 -"M) e - 3c 2l-^l 2 A P ( Tl/3 < M -2 T -(2-a)/ a) 

2 V Vt J \ Vt ) 

= c 3 (l A ^M) (l A ^M) e -3 C2 |-,PA. (2 . 5) 
We now conclude from ((23]), (EHI) and (j2~5l) with 6 = 1 that 

P f,(t, af ,y) = i^lMAf^) > c,t- d ' 2 (l A ^M) (l A *5M) e -3e 2 |,-,PA. 



□ 



The inequality (I2.5P above with 6 = 1/3 will be used later. 
Lemma 2.2 For aZZ M,r,b> 0, there exists Cyi = Ci2(M, r, b) > suc/i that 

P o(^ (0;r) > 6) > C12 > /or a// a G (0, M]. 

Proof. It follows from Lemma 12. II that 

&f Po(r| (0 , r) > 6) = mf / 1^(6,0,^ 

ae(0,M] ; ae(0,M] 7s(0,r) 

> ci-^f-LAll / fM Al ^- Cl 

VV6 JJB(0,r)\Vb J 



e ^ ' b dy. 



□ 



Lemma 2.3 Suppose that M and r are positive constants. Then there is a constant C\ 3 
C 13 (M, a, d, r) £ (0, 1/3] such that for all a € (0, M] and u,v £ R d , 

PB(u,r)UB(v,r) (V 3 , U -> V ) ^ C 13 (J a (u, v) A 1) . 
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Proof. If \u — v\ < r/2, by Lemma 12. II 
PB{u,r)uB { v,r)(V^,v) > ijf ^ (1/3, 0, *) > d (rAlfe"^ > c 3 > <*(./*(«,«) A 1). 



Let E = B(u,r) U B(v,r). If |it - v\ > r/2, with £i = B(u,r/8) and £ 3 = B(v,r/8), we have by 
the strong Markov property and the Levy system (jl.2p of X a that 



p|(l/3, u, «) > E n [p|(l/3 - , , w ) : < 1/3, G £3 

^ Pe 1 { s i u i w )(^J J a (w,z)p%(l/3 - s,z,v)dz S j dw^j ds 



> 



□ 



1/3 

inf J a {w,z)) [ ' F u (t Ei > s) ( [ p a E {l/3-s,z,v)dz)ds 

> F u (t Ei >1/3)( inf J a (w,z)) [ [ p% 3 (l/3-s,z,v)dzds 

\weBi,«eB 3 J Jo Je 3 

= ¥ u (t Ei > 1/3) ( inf - z\)) [ 1/3 F v (te 3 > s)ds 

> h u (r El > 1/3) ( inf j a (| W -z|))p,(rI 3 >l/3). 

Thus by Lemma O 

P a B { u, r) uB M mu,v) > I (Po(r^ (0 , r/8) > 1/3)) 2 mf e£3 r(h - *|) 

> c 4 f(\u-v\) > c 4 (J a (u,v) Al). 

Recall that the function h^,(t,x,y) is defined in (|1.4|) . 

Theorem 2.4 Suppose that M and T are positive constants. There are positive constants C{ = 
Ci(M,Ri,r Q , A , X,T,a), i = 14, 15, such that for all a G (0, M] and (t,x,y) e (0,T] x D x D 

p a D (t,x,y)>C lA h a Cl5 (t,x,y). (2.6) 

Proof. Since t _1 / 2 D satisfies the uniform interior ball condition with radius i?i(T) -1 / 2 for every 
< t < T, there exist 5 = 8{R X ,T) G (0, i?i(T) -1 / 2 ) and L = L(R X ,T) > 1 such that for all t G 
(0, T] and x,y G D, we can choose £ x G (t^^D) n B^^x, L5) and £ y G (t^^D) n B(t~^ 2 y, LS) 
with £(&., 25) n B(£„, 25) = and £(&., 25) U flfo, 25) C r 1/2 £>. 
Let st := t~ l / 2 x and j/ t := ^ 1/2 y. Note that by fl23|) with 6 = 1/3, 

/ P ^ 1/2D (l/3,x t ,u)du > Cl ( S -^-Al) f (5 t - 1/2D (u)Al)e- c ^- u \ 2 du 
JB^ X ,S) \ Vt J Jb(£ X) 6) 

> C1 (^Al) ( -*') ,f |Bfe,i)| 

> C3 (^>4 (, 7> 
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Similarly 



J P^ 1/2D (V3, Vt,u)du > c 3 A l) . (2.8) 

Now we deal with the cases \xt — yt\ > S/8 and \xt — yt\ < 8/8 separately. Recall the definition 
of at from (|23|) . 



Case -Z: Suppose — yt| > 5/8. Note that by the semigroup property and Lemma l2~3l 

rf-i/2 D ( 1 /3, x t , u)p at _ 1/2 n (l/3, it, vV 1 !^ (1/3, u, y t )dudv 

JM y ,5)JB(( x ,5) 

Pt-i/2 nUA x t ,u)p a ' s/2)uB(v s/2) (1/3, u, u )p a l 1/2 (1/3, w, yt)dudu 
B{i v ,S)JB{i x ,S) ' ' " 



> 



> 



>C4 / / p a * /2 Jl/3,x t ,u)(J at (u,v)Al) P a t (l/3,v,y t )dudv 



>c 5 ( inf (J at (ti,«)Al) / / p^ 1/2D (l/3,x t ,u)p a \ /2D (l/3,v,y t )dudv. 



It then follows from (|277|) — ([278]) that 

p^ 1/aD (l 1 x t ,Ift) > c 6 ( inf (J^u,*) A 1)) A (M*) A ^ . (2 . 9) 
Using the fact that 

J at (x t ,y t ) = a a t 1+d / 2 ^( ( i,a)|x-y|-( d+a ) = t 1+d / 2 J a (x, y) (2.10) 

and the assumption \xt — yt\ > 5/8 which implies that \u — v\ < 2(1 + L)5 + \xt — yt\ < (17 + 
lQL)\xt — yt\, we have 

, ^ a ,K nu JJ at (^v)Al)>c 7 (J^(x t ,y t )Al)=c 7 (t 1+d / 2 J a (x,y)Al). (2.11) 
Thus combining (|2.9p and (|2.1ip with (|2.ip . we conclude that for |x t — J/t| > S/8 

P a D (t,x,y) = r d / 2 p^ 1/2D (i,t-y 2 x,ry 2 y) 

> c 8 t-^ (^M A ^ f A A {t w r{x , y) A i) 

= ^^Alj^Al)^,^^). (2.12) 
Case 2: Suppose \x t — yt\ < 5/8. By the semigroup property, 

, , 2D (l,x t ,y t ) > / / p^ 1/2D (l/3,x t ,u)p^_ 1/2D (l/3,u,v)p^_ 1/2D {l/3,v,y t )dudv. 

J B(£y,5) Jb^xjS) 

(2.13) 



p at 
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By with b = 1/3, we have for every (u,v) G B(£ y ,5) x B(£ x ,6), 

P t tt ! 1/2D (l/3, u, u) > c 9 (<W D (u) A 1) (^bW A 1) e -«*ol«-«l a > Cll (<5 A l) 2 . 
Thus by <[277 ]) - ([2^]) and P^3D . 

I«W1,«.,*) > Cl ,(^Al)(^Al) 

> CK (^Al)(^Al)( ( '«''j- (l , S )Al). (2.14) 

Combining (|2.12p and (|2.14p with Lemma 12. 11 we have proved the theorem. □ 



3 Upper bound estimate 



In this section, we will establish upper bound estimate for X a in any open set D (not necessarily 
connected) satisfying a weaker version of the uniform exterior ball condition. 

Suppose that U is a C 1 ' open set with C ' characteristics (R, A). Without loss of generality, we 
can always assume that R < 1 and A > 1. By definition, for every Q G dU, there is a C 1 ' 1 -function 
(t> Q : R*- 1 -> R satisfying Q (O) = 0, V0 Q (O) = (0,...,0), HV^qIU < A, |V0 Q (x) - V0 Q (z)| < 
A|x — z\, and an orthonormal coordinate system CSq y = (y, yd) with origin at Q such that 
B(Q, R)HU = {y = % y d ) G B(0, R) in CS Q : y d > Define 

p Q {x) := x d - 4>q{x), 

where (x, Xd) is the coordinates of x in CSq. Note that for every Q G dU and x G B(Q,R) n ?7, 
we have (1 + A 2 ) _1//2 / oq(x) < 5u(x) < Pq(x). We define for ri,r2 > 

^Q(n,r 2 ) := G U : ri > > 0, |y| < r 2 } . 

We recall the following key estimates from [13^ Lemma 3.5]. 

Lemma 3.1 Suppose R G (0,1], M G (0,oo) and A G [l,oo) are constants, and let ro : = 
fl/(4Vl + A 2 ). T/iere are constants S = 8q(R, M, A, a) G (0,r ), C 16 = C 16 (R, M, A, a) > 
such that for all a G (0, M], A > 1, C 1 ' open set U with characteristics (R,A), Q G dU and 
x G Uq(X^ 1 5q, A _1 ro) with x = 0, 



A^a G f/ < CieAMz) 

C/ Q (A-l4 0i A- 1 r ) 



and 



E, 



^(A-^o.A-iro) 



< dsX-Huix). 



(3.1) 

(3.2) 



We note that 

F x (t% > 1/4) < F x 

< 4E 



^Q(<5o,fo) 



> 1/4 + 



AT". G U and t£ M _ x < 1/4 



r U Q (S ,r Q ) 



+ 



UQ(S ,r ) 



G U 



Thus, by f|3TT ]) - f[3T2p with A 
lemma. 



1 and a simple geometric consideration, we obtain the following 
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Lemma 3.2 Suppose that M > and U is a C 1 ' open set with the characteristics (R,A). There 
exists C17 = Cn(A, R, M, a) > suc/i i/iai /or a// a € (0, M] and x G U, 

F x (t% > 1/4) < C 17 Su(x). 

In particular, we have the following. 

Corollary 3.3 Suppose that M and r\ are positive constants and E := {x G R d : \x — xq\ > n}. 
There exists C\% = C%s(ri,M, a) > independent of xq such that for all a G (0, M] and x G E, 

F x (t e > 1/4) < C 18 5 E (x). 

The proof of the next lemma is similar to that of Lemma 2] , which is a variation of the proof 
of [9l Lemma 2.2]. We give the proof here for the sake of completeness. 

Lemma 3.4 Suppose that Ei,Es,E are open subsets ofM. d with Ei,E^ C E and dist(-Ei, Es) > 0. 
For any n > 1, let E 2 i, i = 1, ... ,n ; be disjoint Borel subsets with U™ =1 -E2,i = E \ {E\ U E3). If 
x G Ei and y G E%, then for all a > and t > 0, 



p%(t,x,y)< y>s(X?a 6%) ( sup p|(s,z,y) + (i A E^rfJ) sup J a (u,z) . 



1=1 



(3.3) 



Proof. Using the strong Markov property, we have 

p% (t, x,y) = E x [p% (t - r% 1 , ,y):r El <t 

n 

= E * [pe (* " r El , ,y):r El <t, X^ G E 2 , 



i=l 



+ E X 



p%(t - r% v X\ ,y) : r El < t,X% G E 3 



--: I + II. 



Clearly 



I < VFjr^a?" £E 2 A[ sup p a E (s,z,y) 
~{ K 1 J \s<t,z£E 2tl 



n 



i=l 

On the other hand, by (jl.2p . 

ft 



< VP, [X% eE 2 A[ sup p a E (s,z,y) 



IJ = / (/ p a El (s,x,u)( [ J a {u,z)p a E {t- s,z,y)dz\du\ds 
Jo \Je! \Je s J J 

< ( sup J a (u,z) ) [ F x (t Ei >s)( [ p a E (t-s,z,y)dz)ds 
\ueEi,z£E a J Jo \Je 3 J 

< [ ¥ x (T Ei >s)ds sup J a {u,z) <(tAE x [r El ]) sup J a {u,z). 

JO u£E 1 ,z£E s u£E 1 ,z£E 3 

This completes the proof of the lemma. □ 
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Theorem 3.5 Suppose that M > is a constant and that D is an open set satisfying a weaker 
version of the uniform exterior ball condition with radius R\ > 0. There exists a positive constant 
C19 = Cig(M, R\) such that for all a G (0, M] and x,y G D, 

p%(l/2,x,y)<C 19 (6 D (x)Al) (e'^ 2 ^ + (f (\x - y\) A 1)) . (3.4) 

Proof. First note that for every xq G M d , {z G M. d : \x — xq\ > R\} is a C 1 ' 1 open set with 
characteristics (R, A) depending only on Ri and d. Let ro and 5q be the positive constants in 
Lemma 1331 for U = {z G R d : \x - x \ > Ri}. 
It follows from Corollary 11.21 that 

p%(l/2,x,y) <p a (l/2,x,y) <c x (e'^ ^ + {j a {\x - y\) M)) , 

so it suffices to prove the theorem for x G D with 5d(x) < 5o/(32). 

Now fix x G D with 5d(x) < 5o/(32) and let Q G dD be such that \x — Q\ = 5d(x). Let 
B Q C D c be the ball with radius R x so that dB Q n dD = {Q} and E := (5q) c . Observe that 
5e{x) = 5 D (x) = \x - Q\. 

When \x — y\ < \J dC 2 V ((5o + ro)/2), we have from Corollary 1 1 . 2 1 that 

p a (l/2,x,y) > c 2 e- C3lx - yl2 > c 4 > and sup p a (l/4, z, y) < c 5 . 

z£R d 

Thus, by the semigroup property and Corollary 13.31 

p a D (l/2,x,y) = [ p a D {l/4,x,z)p a D (l/4,z,y)dz 



1 D 

< sup P a D (l/4,z,y)¥ x (r^>l/4) 

z<=D 

< supp a (l/4,z,y)F x (T%>l/4) 

z&R d 

< cq5 e {x) = c 6 5 D (x) < c 7 6 D (x)p a (l/2,x,y) . (3.5) 



Finally we consider the case that \x-y\ > VdC^V ({5 +r Q )/2) (and 5 D (x) < <5 /(32)). 

There is a C^-function <j) : M^ 1 R satisfying 0(0) = 0, V0(O) = (0,... ,0), ||V0||oo < A, 
|V0(u;) — V0(z)| < A|u> — z|, and an orthonormal coordinate system CS with its origin at Q such 
that 

B(Q, R)r\E = {z= (z, z d ) G 5(0, i?) in CS : z d > 4>{z)} 
and that x has coordinate (0,<5z)(x)) in this CS. Let 

£1 := {z = (z, z d ) in CS : < z d - (p{z) < 5 /8, \z\ < r /8} , 

E 3 := {z £ E : \z - x\ > \x - y\/2} and E 2 := E \ (E 1 U E 3 ). Note that \z - x\ > (5 + r )/4 for 
2 G £3. So, if u G Si and z G £3, then 

|m — z\ > \z — x\ — \x — u\ > \z — x\ — (5q + ro)/8 > -\z — x\ > -\x — y\. (3.6) 



Thus 



sup J a (u,z) < sup J a (u,z) 

u£Ei,z£E 3 (u,z):\u-z\>\\x-y\ 
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< f(\x - y|/4) = {f(\x - y|/4) A j M ((<5 + r )/8)) . (3.7) 

If z G E2, then |z — y| > \x — y\ — \x — z\ > \x — y\/2. We also observe that for every /3 > d/4, 
sup s<1/ / 2 s~ d / 2 e~ l3 / s = 2 d l 2 e~ 2 ^ . By Corollary 11.21 and these observations, 

SUp p a (s,Z,y) < C 3 SUp ^/2 e -|^| 2 /(ft S ) + ^/2 AsJ a (X)!/) j\ 
s<l/2,ze-B 2 s<l/2,zGE 2 ^ ' 

= C^e-\-y\ 2 ^ + ^-f{\x-y\/2) 

= c 3 2 d / 2 e~\ x -y\ 2 ^ + ^ {f(\x - 2/|/2) A j M ((5 + r„)/4)) 
< c 8 (e-l— yl 2 ^ 2 + - 7/1) A 1)) (3.8) 

for some cs > 0. Applying Lemmas 13.11 and 13.41 we obtain, 

P %(l/2,x,y) < c 9 ( e -l— yl 2 /(2C 2 ) + ^ _ y{) A ^ (^(x^ g E) + E^rgj) 

< do ( e -|-^l 2 /(2C 2 ) + _ y|) A ^ 

= do^(x) ( e -l^l 2 /(2C 2 ) + _ yj) A ^ 

Therefore 

j&(1/2, x,y) < P %(l/2,x,y) < c 10 5 D (x) (e^'^^ + ^ - y|) A l)) . 



□ 



Theorem 3.6 Assume that M > is a constant and that D is an open set satisfying a weaker 
version of the uniform exterior ball condition with radius R\ > 0. For every T > 0, there exists a 
positive constant C20 = C2o(T, R\,a, M) such that for all a G [0, M] and (t, x, y) G (0, T] x D x D, 



p a D {t,x,y) <C 20 1 



A *?&\ (l A (t-*/ 2 e-\*-y\ 2 /^ + {t^ 2 A tr(x, y))) . (3.9) 



Proof. Fix T,M > and recall that a* := at( 2 ~ a )/( 2cl ) < MT( 2 - a )/( 2a ). Note that r 1 / 2 !) is an 
open set satisfying a weaker version of the uniform exterior ball condition with radius T -1 / 2 i?i > 
for every t G (0, T]. Thus, by Theorem 13. 5| there exists a positive constant c\ = c\ (T, R\ , a, M) 
such that for all t G (0, T] and a G (0, M], 

p?_ 1/2D (l/2, x, y) < d( e -l*-»l a /(2C a ) + (/t ()a . _ y{) A l))^_ 1/2D ( x) . (3.10) 
Thus by ([23]), rf2~T0j) and IjOty . for every t < T, 

P a D (t/2,x, y ) = r d / 2 ^i 1/2D (i/2,r 1 / 2 x,r 1 / 2 y) 

< ci i^ 2 ( e -l*-y| 2 /(2C 2t ) + ( | x _ j,|/ t V2) A 1} ) St _ i/2D{t -l/2 x) 
= d (r^e-l^l 2 /^*) + (^/ 2 a f J a (x, y))) 
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By symmetry, the above inequality holds with the roles of x and y interchanged. Using the semi- 
group property and Corollary 11.21 (twice), for t < T, 



p a D (t,x,y) = / p a D (t/2,x,z)p a D (t/2,z,y)dz 

< c 3 fotofofo) f p"(2Cit,x,z)p«(2Cit,z,y)dz 

£ Ct SdWSdW ^i/^-yf UK-It) + {t - d /2 A y)) 

This with Corollary 11.21 proves the upper bound (|3.9p by noting that 

(1 A u)(l A v) = min{l, u, v, uv} for u, v > 0. 

□ 

We point out that, in view of Theorem 12.41 the above upper bound estimate (|3.9p is sharp when 
x and y are in the same component of D. However it is not sharp when x and y are in different 
components of D, since in this case when a — > 0, it does not go to zero and thus does not give 
the sharp upper bound for the Dirichlet heat kernel p D (t,x,y) of Brownian motion in D. Next we 
improve the above estimate to get the sharp estimate stated in Theorem 13.91 below. 

For the remainder of this section, we continue assume D is an open set satisfying a weaker 
version of the uniform exterior ball condition with radius Ri > 0. It is easy to see that the distance 
between any two distinct connected components of D is at least R* for some R* > that depends 
only on R\. Without loss of generality, we assume that R* = R\. Observe that for cq > 0, r > vq 
and t > 0, 

t - d /2 e -c r Vt < Cl ^/ 2(t/r y/ 2+ l = C1 J_ < Cir a-2 (3 n) 

where C\ > depends only on cq, tq and and d. This implies that for x and y in different components 
of D, the jumping kernel component tJ x (x, y) dominates the Gaussian component t~ d ^ 2 e~^ x ~ y ^ 2 ^ C2t . 
This fact will be used several times in the rest of this section. 

By Theorem I3.6| we only need to consider the case when x and y are in different components 
of D. Recall that , for any x G D, D{x) denotes the connected component of D that contains x. 

First we give an interior upper bound of p D (t,x,y) when x and y are in different components 
of D. 

Lemma 3.7 Assume that M > is a constant and that D is an open set satisfying a weaker 
version of the uniform exterior ball condition with radius R± > 0. For every T > 0, there exists a 
positive constant C21 = C2i(T, R\,a, M) such that for all a G (0, M], t G (0, T] and x, y in different 
components of D, 

p a D (t,x,y)<C 21 a a t\x-y\- d - a . 



Proof. Using the strong Markov property and (jl.2p . we have for t < T, 
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E, 



' P a D (t-r a Dix) ,X^,y) :t« d{x) <*,X% w £D\D(x) 



t 

\J ' D{x) 
a rt 



< CI 



< Ci 



tM a 
a 



p a D{x) {s,x,u) 

p a (s, x, u) 



D\D(x) 



J a (u, z)p a D (t — s, z, y)dz J du J ds 



o \JD(x) 
a rt 



' D\D(x) 



tM 



(t~ d/2 A (i J M (n, z))) p a (t - s, z, y)dz^j du^j ds 
- jf (j P a (s, x, u) (j (t- d ' 2 A (U M (u, *))) p a (t - s, z, y)dz^j du^j ds. (3.12) 

In the second to the last inequality above, we have used the facts that J M (u,z) < j M (R\) and 
t < T. By Corollary [LSI p a {s,x,u) < c 2 p M \C 2 s,x,u), p a (t - s,z,y) < c 2 p M {C\{t - s),z,y) and 
t~ d / 2 A (tJ M (u,z)) < c 2 p M (t,u, z). Thus, using the semigroup property and Corollary 11,21 from 
(|3,12p we obtain that 



p a D (t,x,y) < c 3 



tM c 



\JR d 



p M (Cls,x,u) 



P 



M (t,u,z)p M (C%(t - s),z,y)dz )du)ds 



a 1 



a rt 



C3 



'p M ((Cl + l)t,x,y)ds 



< c A a a t 



tM a 

a ( t -d/2 e -\x-y\y{C 2 {Cl+l)t) + t -d/2 A 



x-y 



\d+a 



-d—a 



< c^a a t\x — y\ 

In the last inequality above, we have used the fact that \x — y\ > R±. 



□ 



Theorem 3.8 Assume that M > is a constant and that D is an open set satisfying a weaker 
version of the uniform exterior ball condition with radius R% > 0. There exists a positive constant 
C 22 = C 22 (M, Ri) such that for all a £ (0, M] and x,y in different components of D. 

p a D (hoo,y) < C 22 (5 D (x) M)(8 D (y) M) (f(\x - y\) M) . 

Proof. We first claim that 

p%(l/2, x, y) < Cl a Q (5 D (x) A 1) (\x - y\- d - a A l) . (3.13) 

Recall that for every xq £ M. d , {z 6 M. d : \x — xo\ > -Ri/4} is a C 1 ' 1 open set with characteristics 
(R, A) depending only on R\ and d. Let ro and 5q be the positive constants in Lemma 13,11 for 
U = {z £ R d : \x - x \ > Ri/4}. It follows from Lemma [3771 that 

p a D (l/2,x,y) <c x (f(\x-y\) A 1) . 

So it suffices to prove (13. 13ft for x 6 D with 5d(x) < <5o/(32). 

Now fix x £ D with Sjj(x) < 5o/(32) and let Q € dD be such that \x — Q\ = 5n(x). Let 
B Q := B(Q, R x /A) C D c be the ball with radius R x /4 so that dB Q n dD = {Q}. 
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There is a C^-function <j) : R^ 1 R satisfying 0(0) = 0, V0(O) = (0,... ,0), ||V0||oo < A, 
\V(j)(w) — \7cj){z)\ < A\w — z\, and an orthonormal coordinate system CS with its origin at Q such 
that 

B(Q, R) n (Bq) c = {z = (z, z d ) G B(0, R) in CS : z d > 0(1)} 
and that x has coordinate (0,5d(x)) in this CS 1 . Let 

E:=DU(B(Q,R 1 /2)\'Bd, 

E l := {z = (z, z d ) in CS : < z d - <f>(7) < d /8, \z\ < r /8} , 

E 3 :={ze E: \z-x\ > \x-y\/2}, E 2 ,i := (E \ {E\ U E 3 ))0 D(y) and E 2 , 2 := E \ (Ex U E 3 U D{y)). 
Observe that 5g^^ c (x) = 5e(x) = 5d(x) = \x — Q\ and \x — y\> R\ > R > ((<5o + ro)/2). So, by 
dS^D- ([gZD , dist(J5i,£7 3 ) > and 

sup J*(u,z) <c 2 (f(\x-y\) Al). (3.14) 

If z G E 2 i,i = 1,2, then |z — y| > |x — y| — |x — z\ > \x — y\/2. Thus by the same argument as 
the one in (|3.8p . if |x — y| > \ / dC 2 , we have by (|3.1ip 

sup p a (s,z,y)<c 3 (e-\ x -y\ 2 / C2 + (j a (\x-y\)Al)) < c 4 (j M (\x - y\) A l) . (3.15) 

s<i/2,^e£2,i v 7 

If \x — y\ < \/dC 2 , since R\ < \x — y\, we also have 

sup P a {s,z,y) < C 3 sup (s- d / 2 e-^ 2 ^ c ^ + (s- d/2 AsJ a {z,y))) 

s<l/2,zEE 2 ,i s<l/2, ze£ 2 ,i ^ ' 

< C 5 SUp ( s -d/2 e -Hi/{2C 2S ) + ( s -d/2 A s ja (JJl) ^ 

s<l/2 ^ 7 

< c 6 < c 7 j M (|x-y|). (3.16) 
On the other hand, since D{y) C E^, by Lemma 13.71 

sup P a E (s,z,y)<C 21 sup sJ a (|z-y|) < c 8 (j a (\x-y\) Al). (3.17) 

s<l/2,^e-B 2 ,2 s<l/2, z€E 2 ,2 



Furthermore, since dist(£i, D(y)) > Ri/2, by the Levy system (jl.2 
¥ X (X% eE 2jl ) < F X (X% (=D(y)) 



p% 1 (s,x,u) / J a (u,z)dz \ du \ ds 

\JE 1 \JD(y) J J 



a a 



/ P% 1 ( s i x i u ) I J 1 (u,z)dz \ du) ds 

Jo \Je 1 \JD(3i) J J 

< a Q I [ J 1 {z)dz) [ ([ p a E (s,x,u)du) ds 

\J{\z\>R 1 /2} / Jo \Je 1 J 

< c 9 a a ^ x [T a El \. (3.18) 
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Applying Lemmas 13.11 and 13.41 and combining (|3.14p ^ (|3.18p . we obtain, 
P%(l/2,x,y) 

< F X (X% eE 2jl )l sup p a E (s,z,y)) +¥ X (X% eE 2 , 2 )l sup P a E (s,z,y)) 

1 \s<l/2,zEE 2 ,i J 1 \s<l/2,z£E 2 ,2 J 

+E x [t Ei ]I sup J a (u,z)) 

< c 9 a a E x [T Ei ] I sup p a (s,z,y) ) +P X (X% G (B^) c ) I sup p E (s,z,y)) 

\s<l/2,z£E 2 ,i J 1 \s<l/2,zt=E 2>2 J 

+E x [t Ei ]I sup J a (u,z)\ 

< c 10 a Q (fe(x) Al) (|x-yr d - Q Al) . 

Therefore we have proved the claim (|3.13p . In particular, we have 

p a D (l/2,x,y) < c 10 a a (5 D (x) A 1) (Vl^l 2 /(2C 2 ) + (j x _ y \-d-a A ^ (3 1Q) 

By symmetry, the above inequality holds with the roles of x and y interchanged. It follows from 
the semigroup property that 

P a D (hx,y) 

p a D (l/2,x,z)p a D (l/2,z,y)dz 

D 



p a D (l/2,x,z)p a D (l/2,z,y)dz + / p a D (l/2,x,z)p a D (l/2,z,y)dz. 
D(x) JD\D(x) 

By applying Theorem 13.51 and ()3.19|) , and then applying Corollary 11.21 (twice) , we have 

p a D (hx,y) < cna a (5 D (x)M)(5 D (y)M) [ p 1 (2Clx,z)p 1 (2Clz,y)dz 

< Cll a a (5 D (x) A l)(S D (y) A l)p\4Cl x, y) 

< c 12 a a (5 D (x)M){5 D (y)M) (eH*-^/^) + ^ _ y \-d-« A ^ 

< c 13 a a (5 D (x) A l)(S D (y) A l)|x - y\- d ~ a 

< cm (S D (x) A l){5 D {y) A 1) (J a (x,y) A 1) . 

□ 

Recall that h^(t,x,y) is defined in (|1.4p . 

Theorem 3.9 Assume that M > is a constant and that D is an open set satisfying a weaker 
version of the uniform exterior ball condition with radius R± > 0. For every T > 0, there exists a 
positive constant C23 = C 2 s(T, R±,a, M) such that for all a G [0, M] and (t, x, y) G (0, T] x D x D, 

p a D {t,x,y) < C 2 zh a {2C2) - 1 {t,x,y). 
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Proof. Fix T,M > and recall that a t = a t^- a ^^ < MT^~ a ^^ 2a \ Note that t^^D is an 
open set satisfying a weaker version of the uniform exterior ball condition with radius T -1 / 2 i?i > 
for every t £ (0, T]. Thus, by Theorem 13.81 there exists a positive constant c\ = c\(T, R\,a,M) 
such that for all t € (0, T], a E (0, M] and x,y in different components of D, 

p a t l 1/2D (l,x,y) < d (j at (|x - y|) A 1) (5 t _ 1/aj ,(x) A 1) (5 t -^ D {y) A 1) . (3.20) 

Thus by (]2.ip . (I2.10p and (I3.20p . for every i < T, a 6 (0, M] and x, y in different components of D, 

p a D (t,x,y) = t-Wp^^t-Wxtt-Wy) 

< Cl r d ' 2 (>(|x - y\/t^ 2 ) A l) (^v^r 1 / 2 ,) A l) (5 t - inD (r l l 2 y) A l) 




Combining this result with Theorem 13.61 we have proved the theorem. □ 



4 Large time estimates 

In this section we assume that D is a bounded C 1,1 open set in M d which may be disconnected, 
and we give the proof of Theorem 11.31 (ii) and (hi) . 

Proof of Theorem 11.31 (ii) and (iii). Let D be a bounded C 11 open set in M. d with d > 1. For 
each a > 0, the semigroup of X a,D is Hilbert-Schmidt as, by Theorem ll.il 

f p a D (t,x,y) 2 dxdy = [ p a D (2t, x, x)dx < d ((2t)- d / 2 A (a a 2t)- d / a ) \D\ < oo, 

and hence is compact. For a > 0, let {Xt' D ' k = 1, 2 . . . } be the eigenvalues of —(A + a a A a / 2 )\D, 
arranged in increasing order and repeated according to multiplicity, and {4^ D ■ k = 1, 2, . . . } be 
the corresponding eigenfunctions normalized to have unit L -norm on D. Note that : k = 

1,2,...} forms an orthonormal basis of L 2 (D;dx). It is well known that when a > 0, A"' D is 
strictly positive and simple, and that 4>1' D can be chosen to be strictly positive on D. It follows 
from |18[ Theorem 1.1(h)] that the function a i— > A"'' is continuous on (0, M], and lim a ^o+ A"'' D = 
\®' D := mini<j<fc A^'^ , where D\, • • • , are the connected components of D and A^'^ is the first 
Dirichlet eigenvalue of —A\e>.. Hence there is a constant c\ = ci(D,a, M) > 1 so that 

1/ci < \ a { D < a for every a G (0, M]. (4.1) 

Using Sobolev embedding (see [THl Example 5.1]), it can be shown that {4>1' D , a 6 (0,M]} is 
relatively compact in L 2 (D;dx). Hence by |18|. Theorem 1.1(h)] and the fact that each A"' is 
simple for a > 0, a — > 4>1' D is continuous in L 2 (D;dx) in a £ (0,M]. Furthermore, as a — > 0+, 
any weak limit of (j)^ D is a unit non-negative eigenfunction of — A|o with eigenvalue . Note 
that such (p may not be strictly positive everywhere on D when D is disconnected; it is strictly 
positive on least one component Dj of D where A^'^ = \®' D . It follows that there is a constant 
C2 = C2(-D, a, M) > so that 

sup / 5d(x)4>1' D (x)dx < C2- (4-2) 
ae(o,M] Jd 
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Recall that pjj(t,x,y) admits the following eigenfunction expansion 

oo 

p a D (t,x,y) =J2e- tx 2 D <f% D (x)<ft D (y) for t > and x,y G D. 

k=l 

This implies that 

/ 5 D (x)p a D (t,x,y)5 D (y)dxdy = f]e- tx " D ( [ S D {x)4>f (x)dx) . (4.3) 

J DxD Z~y \JD J 

Consequently, we have 



5 D {x)p a D {t,x,y)6 D (y)dxdy <e- tx ^' D I 5 D (x) 2 dx (4.4) 

DxD JD 



for all a > and t > 0. On the other hand, since 

^ D [x) = e^ D [ p a D (l,x,y)<ft D (y)dy, 



D 



by the upper bound estimate in Theorem 11.3( 1) and (|4.1|) - (|4.2p that there is a constant c 3 
cs(D, a, M) > so that for every a G (0, M] and x G D, 



<P{ (x)<c 3 5 D (x) 5 D {y)^i {y)dy < c 2 c 3 5 D (x). 

JD 

It now follows from (|4,3p that for every that for every a G (0, M] and i > 

5 D (x)p D (t,x,y)5 D (y)dxdy>e- tx "' n ( / ^(x)<' D (x)drc ' 

DxD \JD 

> e-' A "' D Uj&^^&fdxj = (c 2 c 3 )- 2 e-' A " D . (4.5) 

It suffices to prove (ii)-(iii) of Theorem 11.31 for T > 3. For t >T and x,y £ D, observe that 

p a D (t,x,y) = p a D (l,x,z)p a D (t-2,z,w)p a D (^w,y)dzdw. (4.6) 

Since D is bounded, we have by the upper bound estimate in Theorem ll.3( i) and (|4.4p that there 
are constants q = Ci(D, a, M) > 0, i = 4, 5 so that for every a G (0, M], t >T and x,y £ D, 

p a D (t,x,y) < c 4 5 D (x)5 D {y) / 5 D (z)p a D {t - 2, z, w)<5 D (w)(fedu; < c 5 fo^/My) 6 "*^- (4.7) 

(ii) Assume first that D is connected. Since D is bounded and connected, we have by the lower 
bound estimate in Theorem [L3ji) and (|4.5p that there are constants Cj = Ci(D, a, M) > 0, i = 6, 7, 
so that for every o G (0, M], t >T and x,y £ D, 

p a D (t,x,y) > c 6 5 D (x)5 D (y) / 5 D (z)p a D (t - 2, 2, -w)5 D (w)dzdu; > c 7 <5 D (x)(5 Z )(y)e"* A ?" D . (4.8) 

This combined with (|4.7p proves Theorem II. 3f ii). 
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(iii) Now let consider the case that D is disconnected. Note that it follows from (ii) that for 
every t > 1, x G D and y G D(x), 



p a D (t,x,y) >p a D(x) (t,x,y) >c 8 e 5 D (x)5 D {y). 



(4.9) 



Moreover, the above inequality, (|4.7p and the two-sided estimate in Theorem ll.3l fi) yield that there 
are a constant eg := cg(D, a, M) > 1 such that for every a G (0, M], t > and x G D, 



a,D(x) 



eg 1 e~ tA i S D (x) < F x [T a D(x) >t)< cge-'-V 5 D {x 



and 



,{T a D >t)<c 9 e- tx ^ 5 D {x). 



(4.10) 



(4.11) 



For t > T, x € D and y £ D\ D(x), we have by the boundedness of .D, (|4.10p and the lower bound 
estimate in Theorem 1 1.3fi) that 



E, 



p2> fai) (s,:c,«) / J a (z,w)p" D (i-s,w,!/)(i«) i Us (4.12) 

\JD(x) \Jd\D{x) I I 



> c iu a"j \^J^^p a D ( x) (s,x,z)dz J IJ^^Diy)^ - s,w,y)dw J ds 



cio a I n x 
lo 



SDW>*)M 7 D(v) >t - S ) dS 



> 



> 



-2 n 



ciocg-a- I e" sA i 6n(x) e-V-'**! 5 D {y)ds 
c u , C f 2 a^e-^ W ^" ,DW )fo( a ;)fc(y). 



(4.13) 
(4.14) 



On the other hand, using (|4.10p - (|4.1ip . we have from (|4.12|) that for t > T, x G D and y £ D\D(x) 
Pr>(^ x ^y) < cna a / pV.Js, x, z)<iz / p£>(t - s,uj, ] ds 

JO \-/D(a;) / \JD\D(x) 

< cii a a jf P*^) > s) Pj,(r£ > i - s)ds 



< c n c 2 9 a a / e" sA i fe(x)e-( t - s ) A i 5 D {y)ds 
Jo 



cnc 9 a a (5z)(x)(5D(2/)e" 



1 e 

10 



ds 



(4.15) 

Jo 

< Cll c 2 9 a a te- tx " D 6 D (x)5 D (y). (4.16) 
Finally, by (|4.7|) and the same argument that leads to (|4.15|) . we have that for t > T, x,y £ D(x), 

pb(t,x,y) 

= p a D (x) x -> v) + ^ to (* - ^d(x) > M , y) ; ^ < * 
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< c 12 5 D (x)5 D (y)e- tx "' D{X) + c 12 a a 5 D (x)5 D (y)e- tx i D [ e "i^"™ ^ ds (4.17) 

Jo 

< c 12 5 D (x)5 D (y)e- tx i D(x) +c 12 a Q t5 D (x)5 D (y)e- tx i D . (4.18) 
Combining this with (|4.7p - (|4.9p . f|4.14j) and (|4.16p completes the proof of Theorem [L3l iii) . □ 



Remark 4.1 In general, when passing from (I4.13D to (I4.14D . from (I4.15D to (14.16ft and from ()4.17p 
to (]4.18p . a factor t is needed in order to have the lower estimate and the upper estimate that is 
uniform in a 6 (0, M]. Note that for D having at least two connected components, \^' D ^ > \°^ D 
for every a > 0. Since D is a bounded C 1 ' open set, it has only finite many connected components 
Di,--- ,-Dfc. According to [181 Theorem 1.1], as a — > 0, A"' D converges to mini<j<fc \® ,D] , where 
\®' Dj is the first Dirichlet eigenvalue of —A\£> j on domain Dj. Let jo be such that A^'^ = 
mini<j<fc \®' Dj . Then for every x G Dj , we have inf ^\^ D ( X ) _ \^' D j = liin \°^ D ^ — \°^ D = o. 

Moreover, if D has two connected components D\ and D 2 that are isometric to each other, then 
by jTHl Theorem 1.1], for x 6 D\ and y 6 D 2 , 

lim A?' DW = A?' Dl = A?' Da = lim A"' D{?/) . 

□ 
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